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Abstract 

The R-function, and its generalization the R-function, are of fundamental importance in the 
fractional calculus. It has been shown that the solution of the fundamental linear fractional 
differential equation may be expressed in terms of these functions. These functions serve as 
generalizations of the exponential function in the solution of fractional differential equations. 
Because of this central role in the fractional calculus, this paper explores various intrarelationships 
of the R-function, which will be useful in further analysis. 

Relationships of the R-function to the common exponential function, e ' , and its fractional 
derivatives are shown. From the relationships developed, some important approximations are 
observed. Further, the inverse relationships of the exponential function, e ' , in terms of the 
R-function are developed. Also, some approximations for the R-function are developed. 

1. Introduction 


The R-function [ 1 ] defined as 




tsr((n + M 


and its generalization the R-function [2], 


R g ,M c ’ 0 s 


„=0 r((n + l )q-v) 


( 1 . 1 ) 


( 1 . 2 ) 


are of fundamental importance in the fractional calculus. In this paper our interest will be 
confined to t > c = 0, q > 0, and v < q . Lorenzo and Hartley [2] have determined a variety of 
relationships associated with the R- function, including those involving relationships with the 
circular and hyperbolic functions as well as other advanced functions. A few more relationships to 
advanced functions are also presented in the Appendix of this paper. It has been shown ([1] and 
elsewhere) that the solution of the fundamental linear fractional differential equation 

a D?x{t)+ax(t)= f{t) (1.3) 


may be expressed in terms of these functions. As in the case of ordinary differential equations 
combinations and convolutions of /?-functions are used to express the solutions of systems of 
fractional differential equations. Because of this central role in the fractional calculus, and since 
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R q 0 (a,0, t ) = F q (a, t) this paper explores various intrarelationships of the -function, which will 
be useful in further analysis and application. 

The general character of the /^-function is shown in figure 1. Figure 1 shows the effect of 
variations of q with v — 0 and a = ± 1. The exponential character of the function is readily 
observed (see, q = 1 ). 

The Laplace transform of the /?-function, starting at t = c = 0 , 

s q - a 

is derived in reference [2]. It is also noted that 

R l o (a,0,r)= /?,o(l,0,a/)= . (1.5) 

»= 0 n«+i) 


These relationships will be useful in the analysis that follows. This is special, because in general 



(1.6) 

The following useful relationship however, is shown 


M10,*) - 1 J , r— r ^ lyfrv -x . 

3% + l)?-v) Sr((n + l)?-v) 

(1.7) 

therefore 


R q ,{ ] 'h0,at)=ar'-'R q Xa\0,t). 

(1.8) 

Alternatively this may be written as 


R qv {a, Oj^a^-^R^XlAa^t). 

(1-9) 


In what follows in this paper, intrarelationships between i?-functions of different arguments will 
be developed. 



Figure 1 a. Effect of qon /? 0 (l,0,f), 
v = 0.0, a = 1 .0 
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Figure lb. Effect of q on R q 0 (- 1,0, f), 
v = 0.0, a = -1.0 
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2. Relationships for R m0 in Terms of R l 0 


This section will develop the relation 

R m .o (l 0 > 0 = / fe.o («>0> * )) = /(*") m = 1A3... 

We consider first the even cases, in particular m = 2. We have then 

oo j 2/2+1 


* 2 - 0 ( 1,0,r ) §r(2n+2) r(2) ' r(4) ' r(6) 


t t 3 t 5 
.( -|- . 


4-... . 


Now since 


and 


/? 1o (1A0=Z-7^ = 1 + ^ + -71 + -71+-- 

sr(«+i) r(2) r(3) r(4) 

/ in x ^ (— , t t 2 / 3 

10 ( ’ 1 )_ S r(»+i) r(2) + r(3) ' r(4) + "' 


R 

it is readily seen by substitution that 

«2,. 0.0 , ») = \ K (i.o, >)- R ,. 0 (- iaO)- 2 

An alternative approach to this problem is available through the Laplace transform, 

^ 1 1 ' 


e -e 


= Sinh(t ) . 


( 2 . 1 ) 


( 2 . 2 ) 


(2.3) 


(2.4) 


(2.5) 


L{R lo (l,0,t)}=-4-=i 

5 — 1 l 


5-1 5+1 


= -j^cis(2Kk/2)\ 

2 k=Q 


s - cis{lK k / 2) J 

where cw(^) = cos(0)+ i sin(^). The inverse transform of this equation, of course, yields the 
equation (2.5) result. 


( 2 . 6 ) 


The m = 4 case is now considered, then 


«4.,(1.0,<)=£: 


t 


4n+3 


r 3 t 1 t" 


■ + - 


+ ■ 


^r(4« + 4) r(4) r(8) r(i2) 

Again by substitution it is readily verified that 

R 4 , 0 (1,0, /) = f {R uo (1,0,/)- R l0 (- 1,0, /) + iR l0 (i,0, /)- iR l0 (- i,0, /)) 


= -^cis(2Kk/4)R lo (cis(27rk/4\0j)=-^e Mn R l0 (e M,2 ,0,t), 


k = 0 


k=0 


— ^-Sinhit j-~-Sin(r) 


(2.7) 

( 2 . 8 ) 

(2.9) 

( 2 . 10 ) 


Examining the above solutions, equations (2.5) and (2.8), it is observed that the values of the 
coefficients and the a parameter of the R { 0 -functions lay on the unit circle of the complex plane 
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for their respective m values. This will be validated for the m = 3 case. Using the Laplace 
transform we have 


(i.o.<)}= -jL-. 

s -1 


Now the roots of 5 3 — 1 = 0 are 


5 = l l/3 = cosl 


r 2 Kk > 

. . (2nk'\ 


(2 Kk\ 


+ 1 sin 

= cis 


^ ^ J 

1 l 3 J 


3 J 


it =0,1,2. 


Thus 


s — 1 

=Hcis(2xkli{ ./* »' 

3 Ho (s-m(27r£/3J^ 


Inverse transforming yields 


1 ^ 

/? 3 0 (l,0, r) = c/s( 2k k / 3)(/?, 0 (cis(2,7Z k/3\0,t)) 

3 jt=o 


This may also be written as 


RxoilO, th\t,cis(27C k/3) 

The above results, equations (2.5), (2.8), and (2.13), are now generalized to give 


(2.11) 


(2.12) 


(2.13) 


(2.14) 


(2.15) 


1 / \ 

/? m0 (lAf)= — y\cis(2x k / m) [R i 0 (cis(2n k/m\0,t)), m - 1,2,3,... , (2-16) 

m k =o 

also 

i m - 1 

^»o(lA f ) = -I« iW ' , "« ril(&S '” ) ', m = 1,2,3,.... (2.17) 

m Ho 

Thus any R m 0 function may be written in terms of R l 0 for m = 1,2,3... . Consideration of the 
principle value k = 0, gives the result 

*..o(1.°.<Lo2 -«>', (2.18) 

m 

this is found to be a useful approximation for t > m. 

The results of equations (2.16) and (2.17) may be generalized to arbitrary positive value for the 
parameter a , this more general result is given by 


R 0 (afij)=—Ya il - m) ' m cis(2nk/m)R l0 (a Um cis{2xk/m),0,t), a> 0, m = 1,2,3 (2.19) 

m S 

and 

R 0 (a,Oj)=—Ya {i - m),m e iM/m e aVmcis{2sk,mh , a> 0, m = 1,2,3,.... (2.20) 

m k=0 
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For negative values of the a parameter the following forms apply 



( 221 ) 

( 2 . 22 ) 


3. Relationships for R Um0 in Terms of R 


1,0 


In this section we seek to express Ry mi o\fi Um ,0,f)in terms of where m = 1,2,3,— • 

The initial interest will be 1 / m = 1 / 2 . Then applying the Laplace transform 


^{^1/2,0 i a 1 )} — 

Inverse transforming gives 


1 _ 1/ 2 . 1/2 „l/2 

1 s +a s 


a 


s v2 - a v2 s V2 +a ul s-a s-a 


jR l/2,o( a ' /2 ’ 0 > r )=^l,I/2( a Af)+a^ 1 ,o(fl,0,f), 


or 


*1/2.0 ( fll/2 ’°’ f K^r 1/2 ^1.0 M> 0+ a R 10 (a,0,t). 

It terms of the exponential function 

r u 2 ,o (a xn $,t)= 0 d)' 2 e at +ae a ' . 

The m = 4 case is considered next. The Laplace transform is applied 

1 


i{« l , < . 0 ( a , ' 4 ,0,r)}= 


1/4 , 1/4 1/4 -1/4 1/4 . -1/4 

s + a s —ia s +ia 


1/4 1/4 1/4 1/4 A/4 - 1/4 1/4 , • 1/4 

s —a s +a s — la s +ia 


3/4 , 1/4 1/2 . _l/2 _ 1 / 4 . 3/4 

5 -ha s -ha s -ha 


s-a 


3_ „(3-*)/4 5 */4 


=I £ 


k = o s — a 

Inverse transforming yields the desired result 


or 


* u «(» l ' 4 .o,«)=t‘> (w| ' 4 *»-.( a .o.O. 

k = 0 


(3.1) 

(3.2) 

(3.3) 

(3.4) 

(3.5) 

(3.6) 

(3.7) 

( 3 - 8 ) 

(3.9) 
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The general results are seen to be 



( 310 ) 

(3.11) 


These results are now validated for the m = 3 case. Then 

2 


k = 0 


Laplace transforming this equation gives 

1 


a 


2/3 


1/3 „l/3 

s —a 


1/3 1/3 

a s 

■+ + 


5-/3 5-/3 


_ 2/3 , 1/3 1/3 , 2/3 

5 5 + /3 


3 (a, 0,0 

(3.12) 

s 2 ' 3 

5 

(3.13) 

s-a 


i 

_2/3 1/3 1/3 ' 

a ) 5 —a 

(3.14) 


The results of equations 3.10 and 3.11 are extended to the case of a negative a parameter by use 
of the following expressions 



(3.15) 

(3.16) 


4. Relationships for the Rational Form /? m/ „ 0 in Terms of R l/p 0 


In this part we wish to relate R m/p o (l,0,0 to R Vpfi {a,0,t), where m and p are positive integers. 
This is a generalization of the preceding parts. We start with the example /f 3/2 (l,0,t), applying the 
Laplace transform gives 

Z'{^3/2,o(hO’0} =— 571 7 • (^'0 

s -1 

Now the roots of s 3n - 1 = Oare s = 1 2/3 = ct.s(27r£y)however because of the periodicity of the 
roots equivalent results are obtained from c/.s(27r k 13). Thus equation 4.1 may be written as 

ife, ,,(1.0.1)}= 

A> . A . A 


s U2 -c 0 s 


1/2 


— c, 


1/2 „ 
s — c. 


(4.2) 
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where c k = cis(2xk/3). The values for this problem are 


c n 


c, =-0.5 + 0.866/, 


c 2 =-0.5-0.866/. 


The A k are determined from partial fraction expansion to be 


A = 


1 


A =- 


l 


li 


-0.5-0.866/ 


A,=- 


1 


\ 


3 -0.5 + 0.866/ 


The A k are recognized as A, = c, / 3. This gives the following for the transform 




1 


s 3/2 -l 


cis{ 0) cis(2Jt/3) cis( 4n/3) 

H — T—rr~ ; a + - 


3(s l/2 — cw(o)) 3(s 1/2 — c/s(2^/3)) 3(s l/2 -c/s(47r/3)) 

The inverse transform then is given as 

1 2 

R 3/2 ,o (lA t )= ~ZjCis( In k / 3)i? 1/2 ,0 (cis(2 K k/3),0,t) . 

3 o 

As in the previous parts these results are generalized to the following 


(4.3) 


(4.4) 


1 


m-1 


R m ipfi^,t) =—^ c isi^Kklm)R {lpfi (cis(^Kklm\Q,t) m = 1,2,3... p = 1,2,3... . (4.5) 


m TZ 


k = 0 


Since e' a = c/s(a)this may be written as 

1 m_1 / \ 

">i.o 


(4.6) 


These results also may be generalized to include a nonunity value for the a parameter. The general 
form is given by 


1 m ~ l ( \ 

R m Ip, o (a,0, 0= “X a (, ' m )y m cis ( 2n k/m)R u p0 (a l,m cis( 2 k k/m\0,t), 

m k = 0 

a> 0, /n = 1,2,3..., p = 1,2,3.., 


and 


m-l 


r> / \ ^ V ’ Jlnklm r> ( Aim Jlnklm /\ A 

R m/p,o (a,0,0=— 2, a « A/ P ,ol fl « ArJ, 

m k-0 

a> 0, m = 1,2,3...,/) = 1,2,3... . 


(4.7) 


(4.8) 


The above equations (4.7) and (4.8) allow any rational based, q = ml p , R m/p0 -function to be 
expressed in terms of its basis R u 0 -functions. 
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The results of equations (4.7) and (4.8) are extended to the case of negative a parameter by the 
following equations 


(4.9) 


(4.10) 

These results, equations (4.7) to (4.10), are the most general (direct) relationships to the basis 
function presented in this paper. The R l/p 0 are seen as basis functions for any R m/pt0 . 


^m/p, o( 

a, 0,0= 

1 

m-1 

2>°- 

m)/m cis((2k + \}r/m)R u Ja Um cis{{2k + l)r/ 

3 

O 



*=0 








a> 0, 

m = 1,2,3..., 

p = 1,2,3.., 

and 







R m /p, o( _ 

-a, 0,0= 

1 

m-I 

I"" - 

m)/ m i((2k+\)rt / m) t* ( A i m i({lk+\)x / ) 





AH 

Jt=0 







a > 0, 

= 1,2,3...,/? 

= 1,2,3.... 


5. Relationships for R u 0 in Terms of R m/ 0 


This section develops the reciprocal relation to that formed in the previous section. This form 
will be useful in developing the inverse relationships, which follow in later sections. Consider the 
case for p = 2, m = 1 , then the Laplace transform is 

1 


ifc 2 .o(a'' ! A/)}= 


1/2 . _ 1/2 1/2 , 1/2 
s +a s +a 


A/1 Ail All . 1/2 

s —a s +a 

A/2 „V/2 

s a 

- + - 


s-a 


s-a s-a 

Inverse Laplace transforming gives 


or 


^/2.o(«' /2 AO=0^ /2 ^,oM/)+« 1/2 /? 1 .oM,0- 

The case p = 2, m = 4, is now considered. The Laplace transform is given by 


L{*./2,o(« 1/4 A0}= 


1 *1/2 , A/A A/2 *1/4 1/2, -1/4 

1 s + <2 5 — /a 5+ /a 


_ 1 / 2 1/4 1/2 , 1/4 1/2 -1/4 1/2 , - 1/4 

s —a s +a s —ta s +ia 


s i,2 +a l ' 4 s l +a U2 s U2 +a^ j,gW*s kn 


s 2 -a 


k= 0 


s - a 


Inverse transforming yields 


or 


R u „tp u \o,t)='£a** y ‘R u , 2 (afl,t), 

k = 0 


(5.1) 

(5.2) 

(5.3) 

(5.4) 

(5.5) 

(5.6) 

(5.7) 

(5.8) 


k = 0 
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The above results are generalized in the following forms 



( 5 . 9 ) 

( 5 . 10 ) 


These results will now be tested on the following case p = 2 y m = 3. Then 




k = 0 


— Cl ^3/2,0 ^ ^3/ 2,1 / 2 /?3/ 2,i (#>0 5 ^)* 

Applying the Laplace transform gives 


1 2/3 1/3 1/2 

1 a a s 5 


1/2 1/3 _ 3/ 2 „ 3/2 ^ 3/2 

^ —a s —a s —a s — a 


(5.11) 

(5.12) 

(5.13) 


„ ■ -1/3 „l/2 , 2/3 

s + a s -ha 


1/2 


. yr ^1/3 ^1/2 , ^2/3 ) 9 
* a + a s -ha J 


1/3 


(5.14) 


„I/2 _l/3 

s —a 


(5.15) 


providing a validation point for the general form (equation 5.9). 

The results of equations (5.9) and (5.10) may also be extended to the case for a negative a 
parameter. These results are given as 



(5.16) 


(5.17) 
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6. Relating R m/p0 to the Exponential Function R l o (b,0,t) = e bt 


Using the results of the previous sections, it is now possible to express any R 0 (<z,0, f) in terms 

of R lo {b,0,t)=e bl , for q = ml p (rational). Two results are required, equations (3.1 1) and (4.8). 
For clarity of discussion we rewrite equation (3.1 1) in the following terms 

K, | A<) 

y=o 

= %b p - l - J 0 d, J ' p e b \ b> 0, p = 1,2,3.... (6.1) 

j=0 


Now this result may be directly substituted into equation (4.8) to give 

e ZA a e R iA a e J 

m k=0 j= o 


m = 1,2,3..., /7 = 1,23—. (6*2) 


This now may be written as 



(6.3) 

(6.4) 


Thus, from equation (6.3) the generalized exponential function R m/pt0 may now be expressed as a 
function of fractional derivatives of the common exponential function 


1 m-\ £-} 
k=0 j = 0 


0 d/' p exp{{a p " n e i2,ckp,m )t), 
a>0,m = 1,2,3..., p ~ 1,2,3... . 


(6.5) 


These results, equations (6.3) to (6.5), contain the results of equations (2.19), (2.20), (3.10), and 
(3.20). 

The case for negative a parameter follows a similar development as above. Equation (4.10) is 
written as 


m - 1 


a > 0, m = 1,2,3..., p = 1,2,3... , 


( 6 . 6 ) 
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where CC = (2k + l)/t / m. Now equation (6.1) may be substituted into equation (6.5) to give 



These results, equations (6.3) to (6.5), and (6.7) to (6.9) are the most general (direct) expressions 
for the /^-function in terms of the common exponential function presented in this paper. 

7. Inverse Relationships — Relationships for R l 0 in Terms of R mk 


In this and the following sections inverse relationships expressing 
terms of various /?-functions will be developed. Consider 

the exponential function in 

4-}= 1 =, v 

s - a [s - a j(s + a) 

(7.1) 

s a 

~ s ^_ a 2 + s 2_ a 2- 

Upon inverse transforming we have 

(7.2) 

= /?, 0 (a ,0, r ) = /? 2 , (a 2 ,0, r )+ a /? 2 0 (a 2 ,0, / ) . 

(7.3) 

In similar fashion for m = 4 , we have 


jLa, T_ 1 _ (j+fl)(j-w)(j+ifl) 

s-a (5-<3)(5 + <3)(5-ta)(5 + ia) ’ 

(7.4) 

s 3 +as 2 +a 2 s + a i 

” 4 4 

s-a 

(7.5) 

3 n 3 ~ k 

Y* a s 

c 4 ' 

(7.6) 
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Inverse transforming this result yields 


e a '=R l0 {a,0,t)=j^a 3 - k o4* * 4 ,o(« 4 At). (7.7) 


k = 0 


k = 0 


These results generalize to 


m - 1 


e a ' =R l0 {aAt)=^a m - l - k R mk (a m ,0,t), a> 0, m = 1,2,3,.., (7.8) 

k = 0 

id 

W-I . X 

«-'=« ltO (fl,O,0=£fl" o^^orAf), a > 0, m = 1,2,3... . (7.9) 


and 


Jt=0 


These results are now used to test the m = 3 case. Then 

2 


=^l.o( a ’ 0 ’0=S a2 ’ < ' ^3,*( fl3 A0- 

Thus the Laplace transform is 


A=0 


•=s 


2 a 2 '*5* 5 2 +a5 + a 2 


j-a ^j 3 -a 3 s 3 -a 3 


i 2 +fl5+a 2 


(j-a)(5 2 + as + fl 2 


’) 5-a 

For negative values of the a parameter the following forms apply 


(7.10) 

(7.11) 

(7.12) 


and 


e- fl ' =^ iO (-fl,O,0=£(-ir W «’ ,M "**-j t ((-l) r, a"A0. a>0 ’ m = 1.2,3,..., (7.13) 

Jfc=0 

«,..(-o.0.<)=S(-ir"«- W „</**„.„((- if «’ At), a>0, m = 1,2,3.... (7.14) 


8. Inverse Relationships — Relationships for R l 0 in Terms of R Um0 

In this section we seek to express /? 10 (l,0,f)in terms of R Um 0 , where m = 1 , 2 , 3 — • The initial 
interest will be 1/ m=l/2. Then 

<* jkn-m ^-1/2 { 0 £.1/2 £l £3/2 

tf./2,o (l A 0 “ S / 2 + 17 2 ) " r(i/ 2) + fj) + r(3 / 2) + r(2) + r(5 / 2) + ’" ^ 

- (-ift k,2 - U2 _ r U2 t° t 1 ' 2 t l t V2 

1/2,0 ’ sr(£/2+i/2) r(i/2) r(i) + r(3/2) r(2) + r(5/2) 


and 


(8.2) 
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(8.3) 


Therefore it is easily seen that 

*1,0 / ) = — (* L/ 2,0 (^, 0“ *1/2,0 ("" * )) * 

We now consider the case II m = 1/4, Then 


*1/4,0 


oo - Jc 74—3/4 

\—L 

Sr(*/4+i/4) 


^-3/4 ^-1/2 ^-1/4 ^0 ^ 1/4 ^i/2 

r(i/4) + r(2/4) + r(3/4) + r(4/4) + r(5/4) + r(6/4) + ’"' 


(8.4) 


Now it may be shown by substitution that 

^l.o( 1 A0=^( /? l/4(lA0-/ ? l/4,o(- 1 A0+ l7? l/4.o( I 'A0-^l/4,o(- I ‘A0)- ( 8 ‘ 5 ) 


As in the previous section this may be generalized as 

R l0 (l,0,t)=—Ycis(2Kk/m)R Um0 (cis(27tk/m),0,t) m = 1,2,3... . 
™ k =o 


( 8 . 6 ) 


Remembering that /? lo (l,0,t) = e ' , this equation (8.6) is recognized as a decomposition (of the 
ra-th order m — 1,2,3“*) of the exponential function. That is, each of these functions is more basic 
than the exponential function in that the exponential function may readily be expressed in terms of 
the “fractional exponential (i.e., R Um 0 (<3,0,/))” in closed form (without differintegrating). 

The results of equation (8.6) may be generalized to arbitrary value for the parameter a , this 
more general result is given by 



(8.7) 


The approach to the solution for a negative a parameter will now be demonstrated for the case 
1/m = 1/4. Applying the Laplace transform gives 


L{R, 0 (-a,0,/)}= L{e- at } 

1 1 

_ „ , „ _ / 1/4 ~ 1/4 t/ 1/4 “ 1/4 V 1/4 1/4 \( 1/4 ~ 1/4 \ 

5 + a (5 — c Q a )\s — c t a j[s —c 2 a )\s — c 2 a ) 

where c k = cis{{2k + l)/r / 4). This may be written as 


( 8 , 8 ) 


A) 


A 




-+- 


. „ 1/4 1/4 1/4 1/4 1/4 ^ 1/4 ' 1/4 ^ _ 1/4 

s + £ s —c 0 a s — c x a s —c 2 a s —c 3 a 


• (8-9) 
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The A k are found by partial fractions to be A k = -c k /(4a 3/4 ). The general form is then validated 


and given as 


R lo (lO-at)=R uo (-a,0,t)=e- a ' 

= —^-a {l ~ m)/m cis(a)R Um0 {a Um cis(a\0,t), a> 0, m = 1,2,3,..., 


1 m-1 , \ 

= — ^-a {x - m) ' m e ia R xhnfi {a Vm e ia At), a> 0, /n = 1,2,3,- , (8.10) 


m t= , 


*=0 


where <X = (2& + l^r / m. 

9. Inverse Relationships — Relationships for e a ' = R l o (a,0,t) in Terms of R m/p 0 

Using the results of the previous sections, it is now possible to express e m = R ] 0 (a,0,t ) in 
terms of any R q0 (b, 0,f). for q = ml p (rational). Two results are required, equations (8.7) and 
(5.10). For clarity of discussion we rewrite equation (5.10) in the following terms 

RuJbA,)=f 1 b''-'- l \d!''R r , p , 0 (b\0,t} f»0, r = 1,2,3..., m = 1,2,3... . (9.l) 

7=0 

Substituting this equation into equation (8.7) gives 

((a""c«(A )(, 0 ,l), 
a>0, m = 1,2,3..., (9.2) 


i m - 1 

R i0 (a,0,t)=e°’=— £a°- v -m 

m k=0 ;= 0 


m k=0 j=o 


ds ( A (r - j))R rlmJlm {a r,m cis( X r ), 0 , r ), 

a > 0, r = 1,2,3..., m = 1,2,3..., 


(9.3) 


where A = 2 7lk! m. 
Alternatively 


«, 0 ( a ,o,,)= £ - =-x 


rri jt=rQ 7=0 



<2 > 0, r = 1,2,3..., m = 1,2,3... . 


(9-4) 


For the case of a negative <2 parameter, we substitute equation (9.1) into equation (8.10) to give 

e m =R lo {l,0,-at)=R uo (-a,0,t) 

1 -Xd5(a)X(a ,, "cM«)r i «7/'-/ f„„ 4 a"-d,(a))', 0 ,r) 

a>0, r = 1,2,3..., m = 1,2,3,... , (9.5) 


; £=o 
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where (X = (2k + \)iz / m . Thus we have 



(9.6) 

(9.7) 


The expressions (9.3), (9.4), (9.6), and (9.7) are the most general expressions for the exponential 
function in terms of the general (rational) /?-function R r/p0 (and its fractional derivatives) 
presented in this paper. 

Tables 1 and 2 summarize the key /^-function relationships developed in this paper in a common 
form. Table 1 presents the relationships for positive b parameter in the left-hand side of 
R u v (b, 0, t ) function, while Table 2 presents the relationships for a negative b parameter. 


Table 1. J?-Function Relationships — Positive Parameter 



Eq. Nos. 

1 / \ 
R m o (a,0, r) = - £ a i '~ m)/m cis(X)R l 0 (a 1/m m(A),0,r) 
m k=0 

(2.19) 

R v ,M ,m ^)=Y j a {m - l - k)hn R Wm (a,0,t) 

k= 0 

(3.10) 

1 / \ 
m k-0 

(4.7) 

JSr=0 

(5.9) 

i m-] p-\ 

Ku P .M&t)= -XE a(W "" )/ffl M(p ~ j)l)R iJlp ((a p,m cis(pX)),0,t) 

m k=0 j= 0 

(6.4) 

w-l . . 

/?, o (a,O,r)=e af =X«'"' ,_t ^(« ra A0 

k-Q 

! 

(7.8) 

R lo (l,0,at) = R uo (a,0,t)=e a ' 

= ±£ a O-)'» c i s (A)R„.. t (a"-cK(A)0,r) 

m *=o 

(8-7) 


(9.3) 


For this table a > 0, m = 1,2,3..., /? = 1,2,3..., r = 1,2,3..., A = 2 Kklm 
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Table 2. fl-Function Relationships — Negative Parameter 

Eq. Nos. 


R m o (- a,0,t) = (a Um cis(a\0,t) 

m t=0 

(2.21) 

k=0 

(3.15) 

m k= o 

(4.9) 

Jt=Q 

(5.16) 

i m-\ p - 1 / v 

m k= 0 j=0 

(6.8) 

R„(- a , 0,t)=e- =g(-l)"- | -* 0 "- 1 -* « <a ((-lJ*fl , \0,r) 

k= 0 

(7.13) 

R l0 (l,0,-fl/) = R l0 (- a,0, t) = e~ m 

— - -2a i, -- ) '-cis(a)R ym , 0 {a' , -cis(a\0,t) 
m k=0 

(8.10) 

R { 0 {1,0 -a t)=R l0 (- a fi,t) = e~ a ' 

1 « 1-1 r - 1 y \ 

= £X a(r " y ‘ m)/m CiS ( a ( r - " C * 5 (« ^0,/) 

m k=0 j = 0 

(9.6) 

For this table a > 0, m = 1,2,3..., p = 1,2,3..., r = 1,2,3..., a — ( 2k + \)k / m, 


10. Approximation of the /^-Function 

Various approximations may be developed for the /^-function. A few such approximations 
will be developed here. As suggested in section 2, the principle value in equation (2.17) provides 
the basis of such an approximation. The result, equation (2.18), is generalized to 

MlA'Me', t > q> 0.5 . (10.1) 

The approximation is shown in figure 2. 

The following is an improved approximation when t < 1 and 0 < q < 1 . 

t > 0, 0<«?<1. (10.2) 

This approximation is shown graphically in Figure 3. 





! 
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R(q,0,1,0,t) and Approximation (see text) R(q,0,l,0,t) and Approximation 



Figure 2. Approximation of ^ 0 (l,0,t) by je' 



Figure 3. Approximation of /? 0 (l*0,f) by Eq. (10.2) 
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The above approximations, equations (10.1) and (10.2) may be extended to include 
R 0 (a,0,/)in the following manner. By a simple replacement of the t variable with at we may 

rewrite equation ( 10 . 1 ) as 

tf,. o (l,0,af) - ±e a ' , a>0,l>q>0, at > 1 (10.3) 

Now using equation (1.9) along with equation (10.3) we infer the approximation 

= a>0, q>0, at >1 (10.4) 

When the a parameter of the /^-function becomes negative, a different set of approximations is 
required. The following approximation works well for l<q<2 

R q0 (- 1,0,0 * e- (2 - ?) '*' cos(f -nil) 1 < <7 < 2, t>n! 2, (10.5) 

This approximation is shown in figure 4. 



Figure 4. Approximation of R q 0 (- 1,0,0 by Eq. (10.5) 


Improved approximations may be determined for particular values of q by optimizing the 
A,B,C, and D constants in the following equation 

R v0 (- 1,0, r) = Ae~ B> cos(c(D t-a 12)) \<q<2 . (10.6) 

When the a parameter in equation (10.5) takes on values other than —1 the following approximation 
works reasonably well for values of a not too large 

a ~ Uq e ^^ ,2Sa ' lq ' cos(a llq t-JZ /2), t> 1, l<q<2. (l0.7) 
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Figure 5 graphically shows this approximation for a = —4. 



1 1.5 2 2.5 3 3.5 4 4.5 5 

{ 


Figure 5. Approximation of Rg,o(~ 4,0, 0 by Eq. (10.7) 

Approximations for values of q>2 with parameter a negative, require positive values for the 
argument of the real exponential in the approximation. For example for 2 < q < 3 the following 
approximation 

R q0 (- 1,0, r) ~ e (° V +1 45 ?- 2 O' cos((- ,2q + 1.4)r -x/2), 2<q<3, t> 1 (10.8) 


is presented in figure 6. 



Figure 6. Approximation of R q 0 {- 1,0, f) by Eq. (10.8) 
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11. Discussion 


This paper has presented a variety of relationships relating various /^-functions. A key result was 
that R q 0 , with q = ml p and positive rational, may be written in terms of basis /^-functions R lt p 0 , 
equations (4.7) to (4.10). Also, reciprocal relationships have been developed expressing R l/p0 in 
terms of R 0 , with q-m! p and positive rational, (equations (5.9), (5.10), (5.16), and (5.17)). 

It was also determined that R q 0 , with q = m/ p and positive rational, may be written in terms 
of fractional derivatives of R l 0 -functions (i.e., exponential functions), equations (6.3) to (6.9). 

Further, the /?, 0 (exponential) functions may in turn be written as a function of basis functions 
R u (equations (7.8), (7.9), (7.13), and (7.14)). These results have allowed very general 
relationships to be written relating R l 0 to ^ 0 and its fractional derivatives, with q = ml p and 
positive rational, (equations (9.3), (9.4), (9.6), and (9.7)). 

It is expected that the results presented here should be analytically very useful since the 
R qv - function is the solution or solution basis of many fractional differential equations. It is also 
observed that all of the above relationships are expressed as finite series, the lengths of which 
depend on m and p. 

Finally, various approximations of /^-functions with both positive and negative arguments have 
been developed. Clearly these approximations only hint at the possibilities, and much more is 
possible. 
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Appendix 


The following relationships of the /?-function with advanced functions are an extension of those 
presented in reference [2]. The expansions used for the defined functions are all taken from 
reference [3]. 


The product of the exponential function and the complementary error function is given by 
ex pM erfc{+ yfx )= f = * 1/2.-1/2 (+ x ) ■ 

*=0 r(l + «/2j 


The error function as given by 

erf(x)=ex p(-x 2 )£ : 


2/r+l 


n = 0 


X(n+3/2) 

The expansion for Dawson’s integral becomes 


= /? uo (-l,0,x 2 K-,/2(lAx 2 ). 




(A-l) 


(A- 2) 


04-3) 


Many distributions may be expressed in terms of exponential of powers of x (see [3], p.260). 
Since /?, 0 (a,0,t)= R l 0 (l,0, at)—e ai these distributions may also be expressed as /?-functions. 
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